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We investigate a new mechanism for the cosmological QCD phase transition: inhomogeneous nu- 
cleation. The primordial temperature fluctuations, measured to be ST/T ~ 10 -5 , are larger than 
the tiny temperature interval, in which bubbles would form in the standard picture of homogeneous 
nucleation. Thus the bubbles nucleate at cold spots. We find the typical distance between bubble 
centers to be a few meters. This exceeds the estimates from homogeneous nucleation by two orders 
of magnitude. The resulting baryon inhomogeneities may affect primordial nucleosynthesis. 
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A separation of cosmic phases during a first-order QCD 
transition III could give rise to inhomogeneous nucleosyn- 
thesis H-|. During a thermal first-order phase transition 
in a homogeneous medium bubbles nucleate due to statis- 
tical fluctuations (homogeneous nucleation) . Their mean 
separation at nucleation introduces a scale for isother- 
mal inhomogeneities in the early Universe, which may 
influence the local neutron-to-proton ratio, providing in- 
homogeneous initial conditions for nucleosynthesis. The 
baryon inhomogeneities may survive until the time of 
neutron freeze-out, if the mean bubble nucleation dis- 
tance, d nuc , exceeds the diffusion length of the proton. 
Comparing those scales at the time of the QCD tran- 
sition, assuming a thermodynamic transition tempera- 
ture T c = 150 MeV, gives d nuc > 2 m {§. The causal 
scale is set by the Hubble distance at the QCD transi- 
tion, dn = c/H ~ 10 km. 

The order of the QCD transition and the values of its 
parameters are still under debate. Nevertheless there are 
indications from lattice QCD calculations |rj-[Io[| . For the 
physical masses of the quarks the order of the transition is 
still unclear Quenched QCD (no dynamical quarks) 
shows a first-order phase transition with a small latent 
heat, compared to the bag model, and a small surface 
tension, compared to dimensional arguments We as- 
sume that the QCD transition is of first order and that 
the values from quenched lattice QCD (scaled appropri- 
ately by the number of degrees of freedom) are typical 
for the physical QCD transition. Based on these values 
and homogeneous bubble nucleation a small supercool- 
ing, A sc = 1—T{/T C 10 -4 , and a tiny bubble nucleation 
distance, d nuc ~ 1 cm, follow [|l]J . The actual nucleation 
temperature is denoted by Tf . 

We argue that the assumption of homogeneous nucle- 
ation is violated in the early Universe by the inevitable 
density perturbations from inflation or from other seeds 
for structure formation. Those fluctuations in density 
and temperature have been measured by COBE JlJ] 
to have an amplitude of ST/T ~ 10~ 5 . The effect 
of the QCD transition on density perturbations [ ^3|Jl4| 



and gravitational waves |lq | has been studied previously, 
while we investigate the effect of the density perturba- 
tions on the QCD phase transition here. We conclude 
that a first order QCD transition induces an inhomo- 
geneity scale of a few meters. In comparison with het- 
erogeneous nucleation via ad hoc dirt Jl6[ , we do not in- 
troduce any new, unknown objects. Our findings might 
have interesting implications for precision measurements 
of primordial abundances |||5] . 

First-order phase transitions normally proceed via nu- 
cleation of bubbles of the new phase. When the temper- 
ature is spatially uniform and no significant impurities 
are present, the mechanism is homogeneous nucleation. 
The probability to nucleate a bubble of the new phase per 
time and volume is approximated by T ~ T 4 exp[— S{T)]. 
The nucleation action S is the free energy difference of 
the system with and without the nucleating bubble, di- 
vided by the temperature. 

Nucleation is a very rapid process, compared with the 
extremely slow cooling of the Universe. The duration of 



the nucleation period, Ai nuc , is found to be 
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dS/dt\ tf 

The time tf is defined as the moment when the fraction 
of space where nucleations still continue equals 1/e. The 
heat flow preceding the deflagration fronts reheats the 
rest of the Universe. We denote by Vhe&t the effective 
speed by which released latent heat propagates in suf- 
ficient amounts to shut down nucleations. In practice, 
v dcf < ^hcat < c s , where v^ef is the velocity of the de- 
flagration front and c s is the sound speed |Q. In the 
unlikely case of detonations Wheat should be replaced by 
the velocity of the phase boundary in all expressions that 
follow. 

The mean distance between nucleation centers, mea- 
sured immediately after the transition completed, is 



(2) 



This nucleation distance sets the spatial scale for baryon 
number inhomogeneities. 
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Lattice simulations |],[l0| imply that in real-world 
QCD the energy density must change very rapidly in 
a narrow temperature interval. This can be seen from 
the microscopic sound speed in the quark phase, c s = 
(dp/de)V\ Lattice QCD indicates that 3c 2 (T c ) = 
0(0. 1) Thus, the cosmological time-temperature re- 
lation is strongly modified already before the nucleations, 
due to 



dT 
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where t H = 1/H = {SM^/Sneq) 1 / 2 with e q being the 
energy density in the quark phase. This behavior of the 
sound speed increases the nucleation distance because of 
the proportionality At nuc oc 1/[3^(T/)] Q. 

In the thin-wall approximation the nucleation action 
has the following explicit expression: 
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for small supercooling. Assuming further that c s docs 
not change very much during supercooling, the following 
relation holds for the supercooling and nucleation scales: 
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Here we denote by A a relative (dimensionless) temper- 
ature interval and by At a dimensionful time interval. 
§ = S(Tf) is the critical nucleation action, S = 0(100). 

Surface tension and latent heat are provided by lattice 
simulations with quenched QCD only, giving the values 
a = 0.015T C 3 , I = 1.4T C 4 §. Scaling the latent heat for 
the physical QCD leads us to take I = 3T C 4 . 

With these values for the latent heat and surface ten- 
sion, the amount of supercooling is A sc = 2.3 x 10~ 4 . 
FromEq. (|) it follows that A nuc = 1.5 xlO~ 6 . Substitut- 
ing 3c 2 = 0.1 into Eq. (|), we find At nuc = 1.5 x 10~% 
for the duration of the nucleation period. The nucle- 
ation distance depends on the unknown velocity i>heat in 
Eq. (|J). With the value 0.1 for t>hcat, the nucleation dis- 
tance rfnuc.hom would have the value 2.9 x 10 _6 (iH- One 
should take these values with caution, due to large un- 
certainties in I and a. As our reference set of parameters, 
we take: A sc = 10" 4 , A nuc = lO" 6 , Ai nuc = 10-H H . 

In the real Universe the local temperature of the ra- 
diation fluid fluctuates. We decompose the local tem- 
perature T(t, x) into the mean temperature T(t) and the 
perturbation 5T(t,n). The temperature contrast is de- 
noted by A = ST/T. On subhorizon scales in the radi- 
ation dominated epoch, each Fourier coefficient A(t, k) 
oscillates with constant amplitude, which we denote by 
Ax(k). Inflation predicts a Gaussian distribution, 
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We find fTJ for the COBE normalized fTJ rms temper- 
ature fluctuation of the radiation fluid (not of cold dark 
matter) A™ s = 1.0 x 10~ 4 for a primordial Harrison- 
Zel'dovich spectrum. The change of the equation of state 
prior to the QCD transition modifies the temperature- 
energy density relation, A = c 2 5e/{e + p). We may ne- 
glect the pressure p near the critical temperature since 
p<£ q at T c . On the other hand the drop of the sound 
speed enhances the amplitude of the density fluctuations 
proportional to cj 1 ^ 2 [Q. Putting all those effects to- 
gether and allowing for a tilt in the power spectrum, the 
COBE normalized rms temperature fluctuation reads 
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where ko = (aiJ)o- For a Harrison-Zel'dovich spectrum 
(n = 1) and 3c 2 = 0.1, we find A™ s » 2 x 10~ 5 . 

A small scale cut-off in the spectrum of primordial 
temperature fluctuations comes from collisional damp- 
ing by neutrinos |2^Jl^]. The interaction rate of neu- 
trinos is ~ GpT 5 . This has to be compared with the 
angular frequency c s /cph of the acoustic oscillations. At 
the QCD transition neutrinos travel freely on scales 
l v « 4 x 10~ 6 dH- Fluctuations below the diffusion scale 
of neutrinos are washed out, 
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In Ref. |1J] the damping scale from collisional damping 
by neutrinos has been calculated to be fcj) h = 10 4 H at 
T = 150 MeV. The estimate (§) is consistent with this 
damping scale. We assume /smooth = 10~ 4 cZh- The com- 
pression timescale for a homogeneous volume ~ ^smooth 
is St = 7rZ smoot h/ c s ~ 10~ 3 £h- Since St ^> At nuc the 
temperature fluctuations are frozen with respect to the 
time scale of nucleations. As long as Z smoo th exceeds 
the Fermi scale homogeneous bubble nucleation applies 
within these small homogeneous volumes. This is a cru- 
cial difference to the scenario of heterogeneous nucleation 
where bubbles nucleate at ad hoc impurities. 
Let us now investigate bubble nucleation in a Universe 
with spatially inhomogeneous temperature distribution. 
Bubble nucleation effectively takes place while the tem- 
perature drops by the tiny amount A nuc . To determine 
the mechanism of nucleation, we compare A nuc with the 
rms temperature fluctuation A™: 

1. If A™ s < A nuc , the probability to nucleate a bubble 
at a given time is homogeneous in space. This is the case 
of homogeneous nucleation. 

2. If A™ s > A nuc , the probability to nucleate a bubble 
at a given time is inhomogeneous in space. We call this 
inhomogeneous nucleation. 

The quenched lattice QCD data and a COBE normal- 
ized flat spectrum lead to the values A nuc ~ 10~ 6 and 
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as deflagrations. They merge within Ai coo i if A nuc < 
(wdcf/whcat)A™ s . This condition should be clearly ful- 
filled for our reference set of parameters. Thus the cold 
spots have fully been transformed into the hadron phase 
while the rest of the Universe still is in the quark phase. 
The latent heat released in a cold spot propagates in all 
directions, which provides the length scale 




t. 



FIG. 1. Sketch of a first-order QCD transition in the in- 
homogeneous Universe. At t\ the first hadronic bubbles (H) 
nucleate at the coldest spots (light gray), while most of the 
Universe remains in the quark phase (Q). At the bubbles 
inside the cold spots have merged and have grown to bubbles 
as large as the temperature fluctuation scale. At £3 the tran- 
sition is almost finished. The last quark droplets are found in 
the hottest spots (dark gray). 



A5? 1S - 10 



5 . We conclude that the cosmological QCD 
transition may proceed via inhomogeneous nucleation. A 
sketch of inhomogeneous nucleation is shown in Fig. [l]. 
The basic idea is that temperature inhomogeneities de- 
termine the location of bubble nucleation. Bubbles nu- 
cleate first in the cold regions. 

The temperature change at a given point is governed 
by the Hubble expansion and by the temperature fluctua- 
tions. For the fastest changing fluctuations, with angular 
frequency c s /Z smoot h, we find 
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The Hubble expansion is the dominant contribution, as 
typical values are ic 2 s =0.1 from quenched lattice QCD 
and Atp^tn/St » 0.01 from the discussion above. This 
means that the local temperature does never increase, 
except by the released latent heat during bubble growth. 

To gain some insight in the physics of inhomogeneous 
nucleation, let us first inspect a simplified case. We 
have some randomly distributed cold spheres of diam- 



eter / 
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with equal and uniform temperature, which 



is by the amount A™ S T C smaller than the again uniform 
temperature in the rest of the Universe. When the tem- 
perature in the cold spots has dropped to Tf, homoge- 
neous nucleation takes place in them. Due to the Hubble 
expansion the rest of the Universe would need the time 
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to cool down to Tf. Inside each cold spot there is a 
large number of tiny hadron bubbles, assumed to grow 
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If the typical distance from the boundary of a cold spot 
to the boundary of a neighboring cold spot is less than 
^hcat, then no hadronic bubbles can nucleate in the in- 
tervening space. In this case the nucleation process is 
totally dominated by the cold spots, and the average dis- 
tance between their centers gives the spatial scale for 
the resulting inhomogeneities. In the following analysis 
for a more realistic scenario we concentrate in this case, 

^hcat ^ ^smooth- 

The real Universe consists of smooth patches of typ- 
ical linear size Z sm0 oth, their temperatures given by the 
distribution (||). As discussed above, the merging of tiny 
bubbles within a cold spot can here be treated as an in- 
stantaneous process. The fraction of space that is not re- 
heated by the released latent heat (and not transformed 
to hadron phase), is given at time t by 



/(t)«i- / r lhn (t')v(t,t')dt', 



(12) 



where we neglect overlap and merging of heat fronts. At 
time t heat, coming from a cold spot which was trans- 
formed into hadron phase at time t' , occupies the volume 
U(M') - (47r/3)[Z smo . oth /2 + Wheat (t - t')} 3 . The other 
factor in Eq. (|l2|), r^n, is the volume fraction converted 
into the new phase, per physical time and volume as a 
function of the mean temperature T — T(i). is pro- 
portional to the fraction of space for which temperature 
is in the interval [Tf, 1/(1 + dA)]. This fraction of space 
is given by Eq. (§) with A = Tf/T — 1. Rewriting dA by 
means of Eq. (0) leads to the expression 
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where the relevant physical volume is V smoo th = 
(4^/3)(Uooth/2) 3 . 

The end of the nucleation period, tihn, is defined 
through the condition J(tii m ) = 0. We introduce the 
variables N = (1 - T f /T)/A™ s and = N(t ihn ). Since 
c s may be assumed to be constant during the tiny tem- 
perature interval where nucleations actually take place, 
we find from Eq. §): 1 - t/t ihn « 2/(3 C 2)A^ ms (A^ - AO- 
Putting everything together we determine J\f from 
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The COBE normalized spectrum gives 'heat /'smooth = 

2Whoat(3 C 2)- 1 /4( fc/fco )(«-l)/2^ p or Z hoat / Uooth = 

1, 2, 5, 10 we find TV w 0.8, 1.4, 2.1, 2.6, respectively. 

The effective nucleation distance in inhomogeneous nu- 
cleation is defined from the number density of those cold 
spots that acted as nucleation centers, rf n uc,ihn = n^ 1 ^ 3 . 
We find 
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T ihn {t)dt 
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= {-[l-cri(JV/V2)}}~ 1 / 3 I 



smooth- 



(15) 
(16) 



most probable range of cosmological and QCD parame- 
ters. 
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With the above values 'heat /'smooth 

rfnucihn = 1-4, 1.8, 3.0, 4.8 X 'smooth 



1,2,5,10 we get 
where 'smooth ~ 1 m. 



For a COBE normalized spectrum without any tilt and 
with a tilt of n — 1 = 0.2 (where (fc smoo th/fco) 01 ~ 25), 
together with 3c;? = 0.1 and Wheat = 0.1, we find the esti- 
mate 'heat/'smooth ~ 0.4 and 9, correspondingly. Notice 
that the values of Wheat and 3c 2 s are in principle unknown. 
Anyway, we can conclude that the case 'heat > 'smooth is 
a realistic possibility. 

With 2« hoat (3 C 2)- 1 /4( 10 -4 dH / /smooth) < i an d with- 
out positive tilt we are in the region Zheat < 'smooth, 
where the geometry is more complicated and the above 
quantitative analysis does not apply. In this situation 
nucleations take place in the most common cold spots 
(Af ~ 1), which are very close to each other. We ex- 
pect a structure of interconnected baryon-depleted and 
baryon-enriched layers with typical surface 's moo th an< ^ 
thickness Zdcf = UdefAt coo i- In between cZ n uc,hom would 
be the relevant length scale of inhomogeneities. An accu- 
rate analysis of this case requires computer simulations, 
which is beyond the scope of the present work. However, 
it is clear that the result will be different compared with 
homogeneous nucleation. 

We emphasize that inhomogeneous and heterogeneous 
nucleation pi are genuinely different mechanisms, al- 
though they give the same typical scale of a few meters by 
chance. If latent heat and surface tension of QCD would 
turn out to reduce A sc to, e.g., 10~ 6 , instead of 10~ 4 , the 
maximal heterogeneous nucleation distance would fall to 
the centimeter scale, whereas on the distance in inhomo- 
geneous nucleation this would have no effect. 

We have shown that inhomogeneous nucleation dur- 
ing the QCD transition can give rise to an inhomogene- 
ity scale exceeding the proton diffusion scale. The re- 
sulting baryon inhomogeneities could provide inhomoge- 
neous initial conditions for nucleosynthesis. Observable 
deviations from the element abundances predicted by ho- 
mogeneous nucleosynthesis seem to be possible in that 
case §,§. 

In conclusion, we found that inhomogeneous nucleation 
leads to nucleation distances that exceed by two orders of 
magnitude estimates based on homogeneous nucleation. 
We emphasize that this new effect appears for the (today) 
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